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cluster states 

Chol-Ung Choe^'^, Thomas Dahms\ Phihpp H6vel\ and Eckehard Scholl^B 

^ Institut fiir Theoretische Physik, Technische Universitat Berlin, 10623 Berlin, Germany 
^Department of Physics, University of Science, Unjong- District, Pyongyang, DPR Korea 

We study synchronization in delay-coupled oscillator networks, using a master stability function 
approach. Within a generic model of Stuart-Landau oscillators (normal form of super- or subcritical 
Hopf bifurcation) we derive analytical stability conditions and demonstrate that by tuning the 
coupling phase one can easily control the stability of synchronous periodic states. We propose 
the coupling phase as a crucial control parameter to switch between in-phase synchronization or 
desynchronization for general network topologies, or between in-phase, cluster, or splay states in 
unidirectional rings. Our results are robust even for slightly nonidentical elements of the network. 
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Over the last decade, control of dynamical systems and 
stabilization of unstable states have become a central is- 
sue in nonhnear science [H. In parallel, the study of 
coupled systems ranging from a few elements to large 
networks has evolved into a rapidly expanding field [2|. 
To determine the stability of synchronized oscillations 
in networks, Pecora and Carroll introduced a technique 
called master stability function (MSF) which allows 
one to separate the local dynamics of the individual nodes 
from the network topology. Although some recent ap- 
proaches have tried to extend this theory in the presence 
of time delays 0, [B| , up to now control and design of dy- 
namic behavior in complex networks with time delay is 
still in its infancy. 

In this Letter, we aim to fill this gap by developing 
analytical tools for a large class of delay-coupled net- 
works and deriving analytical conditions for controlling 
the different states of synchrony. We identify the cou- 
pling phase as a crucial control parameter and demon- 
strate that by adjusting this phase one can deliberately 
switch between different synchronous oscillatory states of 
the network. We use a generic model describing a wide 
range of systems near a Hopf bifurcation, which allows 
for an analytical treatment, including the calculation of 
the Floquet exponents. These results promise broad ap- 
plicability, since the presence of time delays is of crucial 
importance in a variety of physical, biological, technolog- 
ical, social, ecological, or economic networks where they 
occur, e.g., as propagation delays in communication net- 
works and laser arrays B, 0, electronic circuits 
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We consider A^-dimensional networks of delay-coupled 
Stuart-Landau oscillators {j — 1, . . . , N) 



N 

n=l 



(z„(t-r)-z,(<)), (1) 



the odd-number limitation of time-delayed feedback con- 
trol [l3| and in anticipating chaos synchronization [l^ . 
The topology of the network is determined by the real- 
valued adjacency matrix A = (a^n). Nonzero diagonal 
elements, for instance, correspond to networks with de- 
layed self- feedback. In the following, we consider only 
constant row sum ^ = J^n '^j" ^^^^ that each node is 
subject to the same input for complete synchronization. 
This generalizes the common assumption of zero row sum 
in the MSF approach. The local dynamics of each ele- 
ment is given by the normal form of a supercritical (— ) 
or subcritical (+) Hopf bifurcation: 



/ i^j) = [A + icj =F (1 + il)\Z] 



(2) 



with Zj = rjC^'^^ £ C, time delay t, and complex cou- 
pling strength cr = /ve*'^. Such phase-dependent cou- 
plings have been shown to be important in overcoming 



with real constants A, a; ^ 0, and 7. This system arises 
naturally as a generic expansion near a Hopf bifurcation, 
and is therefore often used as a paradigm for oscillators. 

In the following, we focus on synchronous in-phase, 
cluster, and splay states with a common amplitude rj = 
'ro,m, and phases given by (pj = flmt + jA4>m with A(j>m = 
27rm/7V. The integer m determines the specific state: in- 
phase oscillations correspond to to = 0, while cluster 
and splay states correspond toTO=l,...,A^ — 1. The 
cluster number dc, which determines how many clusters 
of oscillators exist, is given by the least common multiple 
of TO and N divided by m. dc = N corresponds to a 
splay state (l9l |. With the above notation, we obtain for 
in-phase oscillation in general networks and for splay and 
cluster states in ring configurations 

/ N \ 

rQ „^ = ± I A - fj.K cos /3 + if ^ ajn cos j (3a) 

N 

il„i ^ UJT 7^o,m ~ f-^I^ sin + X ^ Ujn sin <i>„,m (3b) 

n=l 

as invariant solutions of ro,m and (pj using the abbrevia- 
tion $„.m = P — + {n — j)A(/),n, which is indepen- 
dent of J in the cases mentioned above The following 
discussion focuses on the supercritical case (upper signs), 
but a similar argument holds also for the subcritical Hopf 
normal form (see discussion at the end of this Letter) . 
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FIG. 1: (Color online) Collective frequency Slo (left) and 
squared amplitude Tq (right) of in-phase oscillation (m = 0) vs 
time delay r for different amplitude of the feedback strength 
(/iA' = 0.3 and 0.7 in (a),(b) and (c),(d), respectively). Black 
(blue) and gray (yellow) curves correspond to a feedback 
phase /3 = and P = flor, respectively. Unphysical solu- 
tions {tq < 0) are dashed. For f3 = Qqt the curves in (a), 
(c) have the same shape, but no unphysical solutions occur. 
Parameters: A = 0.1, a; = 1, 7 = 0. 

Figure [T] shows solutions of Tq (ro = ro,o) and fio for 
in-phase oscillations (to = 0) according to Eqs. ([3]) in de- 
pendence on the time delay r for fixed feedback strength 
liK = 0.3 and 0.7 in panels (a),(b) and (c),(d), respec- 
tively. The black (blue) lines show the behavior for the 
coupHng phase 13 = 0. The collective frequency fio is 
distributed around the intrinsic frequency uj = 1, where 
multiple solutions are obtained with increasing time de- 
lay r. This behavior becomes more pronounced for higher 
fiK (c). The collective amplitude also shows multival- 
ued behavior; spurious solutions with r§ < 0, which 
correspond to amplitude death, are indicated as dashed 
curves. For a coupling phase f3 = floT, these unphysical 
solutions do not occur since > A = 0.1 as shown by the 
gray (yellow) curves in (b),(d). Note that for /3 = Qqt 
the shape of the flo curve in (a),(c) is unchanged, but 
now all points are valid solutions. 

Considering small deviations Srj and S(pj, i.e, rj = 

yields a variational equation for the synchronized state 

i = lN&) [Jim - K^m)C + K{A ® UnMt - t) (4) 

with the 2iV-dimensional vector ^ = (^i, . . . , ^at)"^, 
the N X N identity matrix I^v, and matrices = 

V J2n "-jn sin J2n ""jn COS ) ' ^ 

/ cos - sin \ -fT ^ ( T2rg „ \ 

i^sin$„,„ cos$„,™ y ' 0.™ \T^^rl^{))' 
which is an important generalization of the usual MSF 
approach. 

In order to derive an analytical expression for stability, 
Eq. ([4]) has to be diagonalized in terms of A. To succeed, 
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FIG. 2: (Color online) (a), (b) Master stability function in 
the (Rei^, Im:/)-plane for m = 0, /3 = 0. The gray scale 
(color code) corresponds to the largest real part of the Floquet 
exponents for a given value of the product Kv. All eigenvalues 
of the coupling matrix A for a unidirectional ring lie on the 
black circle. Parameters: (a) = 0.3, r = 2-k, (b) = 
0.08, r = 0.527r, others as in Fig[Tl (c) Distribution of relative 
phases 6 = (fi — 6 around the order parameter for 200 slightly 
nonidentidal elements with different standard deviations cr^ 
of the frequencies u for fi = Qot. Dotted (black) curve: P = 
iloT + n (desynchronization). Other parameters as in (b). 

the rotational matrix Rn,m must not depend on n. This 
is achieved in two cases: (i) By considering only in-phase 
synchronization (m = 0), or (ii) by considering special 
network configurations. 

In case (i), the matrix R„^o = R- with $„^o = ^"0 does 
not depend on n and with J^^ = Eq. (0]) simplifies to 

i^lN(^iJ^ - ^iKR)^ + K{A ® R)^{t - t) (5) 

Diagonalizing A, we arrive at the block-diagonalized vari- 
ational equation: 

Ut) = Jja(i) - KIi[fiCkit) ~ VkCkit ~ r)] , (6) 

where is an eigenvalue of A, fc = 0, 1, 2, . . . , iV— 1, and 
= IX corresponds to the dynamics in the synchroniza- 
tion manifold. Since the coefficient matrices in Eq. ^ do 
not depend on time, the Floquet exponents of the syn- 
chronized periodic state are given by the eigenvalues A 
of the characteristic equation 

det {JJ - AI2 + X ( - Ai + Vke-^^)K] = 0. (7) 

Figures [2ta),(b) depict the MSF, i.e., the largest real 
part of the Floquet exponents, calculated from Eq. ([7]) 
for different coupling parameters. Note that for a uni- 
directionally coupled ring all eigenvalues, i.e., = 
exp(27ri/c/7V) with fc = 0, 1, . . . , iV — 1, are located on the 
black circle. Hence, for the choice of parameters in panel 
(a) all eigenvalues lie in the region of negative maximum 
real part of the Floquet exponent (stable in-phase solu- 
tion), whereas the parameters in panel (b) do not allow 



for synchronization of the unidirectional ring. Further- 
more, it can be shown using Gerschgorin's disk theorem 
[l^ that the eigenvalues are located on or inside this 
circle 5'(0,//) centered at with radius fj, for any net- 
work topology without self-feedback (diagonal elements 
ttjj ~ 0) . The same holds for the circle S {ajj , fi) cen- 
tered at Qjj if self-feedback (constant ajj , j = 1 . . . , N) 
is added while keeping the constant row sum condition. 
Note that the MSF is symmetric with respect to a change 
of sign of Im(z^). 

For coupling phases /3 = floT + 21t: with integer Z, i.e., 
<i>o = 2Z7r, the characteristic equation ([7]) for the Floquet 
exponents A factorizes: 

= {Kvkc-^^ - Kfi T 2r^ - A) (AVfee""^^ - A> - A) . 

(8) 

In the supercritical case (upper sign), the dominant Flo- 
quet exponent is determined by the second factor in 
Eq. ([HI) which gives Re A < for K > for any eigenvalue 
Vk on or inside 5(0, /i), taking into account Tq > 0, and 
hence stable in-phase synchronization for any network 
topology without or with self-feedback, observing the 
constant row sum conditions. A similar equation arises 
for /3 = rioT + {21 + l)7r with fio obtained from Eq. l(3bl) . 
and it can be shown by analogous arguments that there 
exist exponents A with positive real part, which results 
in desynchronization. In conclusion, the synchronous 
(in-phase) dynamics can be stabilized or destabilized by 
proper choice of the coupling phase /3. These results are 
robust even if slightly nonindentical elements are consid- 
ered. Fig. Wic) shows numerical simulations for all-to-all 
coupling of 200 elements with Gaussian frequency distri- 
butions around u — 1 with different standard deviations 
(Ju,- For /? = r^oT the relative phases of the individual os- 
cillators around the order parameter i?e*® = 
(in rotating coordinates) are distributed around a max- 
imum (in-phase) with different sharpness according to 
cr^j, whereas for (3 = -\- tt the phase is uniformly dis- 
tributed regardless of (desynchronization) . 

We now consider case (ii), i.e., special network configu- 
rations, and exemplarily choose a unidirectional ring with 
0'i,i+i = ajv,i = l= M ciiid all other a.ij = 0. Then the 
matrix R„,,„ = R™ with $„^,„ = $,„ = /? - QmT + A4>m 
does not depend on n and the diagonalization of Eq. ([4]) 
yields the same form as Eq. ([6]) with /i = 1, and R re- 
placed by Rm. The eigenvalues of A are explicitly given 
by vk = e^''^'^/^, fc = 0, 1, . . . ,iV - 1, and the Floquet 
exponents A can be calculated from the corresponding 
characteristic equation 

Figure [Sja) shows the stability boundaries of different 
dynamical scenarios in the {K, T)-plane for unidirectional 
coupHng of -/V = 4 oscillators. The coupHng phase is fixed 
at /3 = 0. The gray scale (color code) indicates regions 
of different multistability of in-phase (to = 0), 2-cluster 
(m = 2), and splay states (to = 1, to = 3): Black (blue), 
dark gray (red), light gray (green), and yellow (white) 
color corresponds to regions where one, two, three, or 
four of these dynamical states are stable, respectively. 
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FIG. 3: (Color online) (a) Stability diagram of unidirectional 
ring of Ai' = 4 oscillators in the (/T, T)-plane (/3 = 0). Solid, 
dash-dotted, dashed, and dotted boundaries correspond to 
a stability change of in-phase (?n = 0), 2-cluster (m — 2), 
splay states with m — 1, and m — 3, respectively. The black 
(blue), dark gray (red), light gray (green), and white (yellow) 
regions denote multistability of one, two, three, and four of the 
above states, (b)-(e) Time series of the phase differences for 
a unidirectional ring of four slightly nonidentical oscillators: 
(b) (3 = nor, (c) ^ir - 7r/2, (d) ^ar - tt, (e) Q3T - 3n/2, 
with Qo = 1, f!i = 0.83903, fia = 1, and ^3 = 1.16097. Black 
(blue), dark gray (red), and light gray (green) lines denote 
the differences ip2 ~ fi, ips — ifii, and ip4 — fi, respectively 
(in (b),(d) black (blue) is hidden behind light gray (green)). 
Parameters: A = 0.1, 7 = 0, if = 3, r = 37r, a;i = 0.99757, 
UJ2 = 0.99098, UJ3 = 1.01518 and UJ4 = 0.99496. 



Let us now consider the effects of the coupling phase 
/?. The specific choice of /? = Hqt, — 7r/2, — tt, 
and flsT — 3tt/2 enlarges the stability regime of the 
in-phase, splay (to = 1), cluster, and splay (m ~ 3) 
states, respectively, to the complete (K, T)-plane. This 
can be understood as follows. For $m = 2Z7r, i.e., 
P = flmT — A(j)„i + 21t: with integer Z, the character- 
istic equation can again be factorized as Eq.®. For 
the supercritical case, taking into account Tq > at 
/? = O.m''' — + 21t:, it follows again that the domi- 
nant Floquet exponents have negative real part for any 
K > Q and r. Therefore the unidirectional ring con- 
figuration of Stuart-Landau oscillators exhibits in-phase 
synchrony, splay state and clustering according to the 
choice of the control parameter (i = fioT, fiir — 2tt/N , 
or 0.mT — 2'Kraj'N {ra > 1,7V > 2), respectively, for any 
values of the coupling strength and time-delay. 

To illustrate this further and demonstrate the robust- 
ness of our stability results for slightly nonidentical oscil- 
lators, we choose a set of control parameters AT = 3 and 
T = 37r, denoted by the black cross in Fig-jS^a), for which 
multistability of all four possible synchronization states 
is found for the coupling phase = 0. Figures [3l[b)- 
(e) show time series from numerical simulations of four 
Stuart-Landau oscillators in a unidirectional ring config- 
uration with slightly different frequencies w. For each 
choice of /? in panel (b) - (e) the solutions were ob- 
tained by solving Eqs. ^ such that the solution of 
closest to unity was chosen. The differences of the phases 
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FIG. 4: (Color online) Real part of the Floquet exponents A in 
the subcritical case vs reduced coupling strength (A*' — l)K in 
the all-to-all and star coupling configurations for different A*'. 
The solid (blue) lines denotes ReA inside the synchronization 
manifold for both all-to-all and star coupling. The dashed 
(red) lines show the largest transversal ReA for different A'' 
and all-to-all coupling. Floquet exponents for star coupling 
are independent of A''. Parameters (as in [i3|): A = —0.005, 
cj = 1, 7 = -10, r = 27r/(l - 7A), /3 = 7r/4. 

(fi {i = 2, 3, 4) relative to the first oscillator phase tpi are 
plotted. After transients (note that the transient oscilla- 
tions are not resolved on the time scale chosen), the os- 
cillators behave exactly as predicted by our theory, i.e., 
they lock into in-phase synchronization for [3 = lot (b), 
into a splay state for /3 = lot — n /2 (c), into a 2-cluster 
state for (3 = lot — tt, where ipi = Lp^ and Lp2 = (d), 
and again into a splay state, albeit with inverted ordering 
of the phases, for (3 = lot — 37r/2 (e). 

Finally, for the subcritical Hopf normal form, it can be 
shown that the periodic orbit, which is unstable in the 
uncoupled case, can be stabiHzed in-phase synchronously 
by, e.g., bidirectional ring, star, or all-to-all coupling 
without self- feedback. In these cases, both the synchro- 
nization manifold and the transversal modes are stable. 
There, Floquet exponents A satisfy again Eq. ^ with 
vq = /i. For proper coupling phases /?, we find a finite 
interval of feedback gain K for which the real parts of all 
Floquet exponents are negative. Note that the in-phase 
synchronization manifold coincides with a single oscilla- 
tor with delayed self-feedback as considered in Ref. 
to refute the alleged odd number Hmitation. 

For all-to-all coupHng the adjacency matrix is given 
by ai^i = and a^.j = 1 for i ^ j and i,j = 1, . . . ,N , 
while for star coupling all j = 0, except ai^i = 1 and 
Oj^i = iV - 1 for i = 2, . . . , A^. Figure H shows Re A 
in the subcritical case as a function of coupHng strength 
nK for all-to-all and star coupling. The soHd (blue) lines 
corresponds to Re A inside the synchronization manifold 
(vq = fi = N — 1 for both coupling configurations). For 
all-to-all coupling the transversal eigenvalues of A are 



given by t'fc = — 1 for fc = 1, . . . , A''—! and the correspond- 
ing largest Re A are denoted by the dashed (red) lines in 
Fig. m for different N. For star coupHng, the transversal 
eigenvalues of A are ft = for k = 1, . . . , N — 2 and 
vn-i = —{N — 1), and the corresponding largest Re A 
are marked schematically. We stress that for both all-to- 
all and star coupling there exists an interval of feedback 
strength K in which all Re A < 0. Thus, time delayed 
coupHng results in stabiHzation and in-phase synchro- 
nization. 

In conclusion, we have shown that by tuning the cou- 
pling phase in delay-coupled networks one can easily con- 
trol the stability of synchronous periodic states, and we 
have specified analytic conditions. In general networks 
in-phase synchronization or desynchronization can be 
chosen, and in unidirectional rings either in-phase, clus- 
ter or splay states can be selected. The coupling phase 
is a parameter which is readily accessible, e.g., in optical 
experiments [H. Our results are robust even for slightly 
nonidentical elements of the network. 
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